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SUMMARY 



A generalized ' method of analyzing experimental obser- 
vations in problems of elastic stability is presented in 
which the initial readings of load and deflection may be 
taken at any load less than the critical load. This anal- 
ysis is an extension of a method published by Southwell 
in 1932, in which it was assumed that the initial readings 
are taken at zero load. 



INTRODUCTION 



In reference 1, Southwell presented a method for the 
analysis of experimental observations in problems of elas- 
tic stability. Briefly, the method is concerned with the 
interpretation of simultaneous readings of load and deflec- 
tion. As therein presented, the method requires that the 
initial deflection reading be taken at zero load. In the 
vicinity of zero load, deflection readings are usually 
somewhat questionable. The deflection readings are reli- 
able only after enough load has been applied thoroughly to 
seat the specimen and the loading fixtures. Furthermore, 
it is not always convenient to take the initial deflection 
reading at zero load. Something may also happen to render 
the first few readings valueless and it may not be practi- 
cable to repeat them. For use under such circumstances* 
a more general method has been devised wherein the initial 
readings may be taken at any load less than the critical 
load. 

The general method of analyzing experimental observa- 
tions in problems of elastic stability is presented in this 
paper. Reference 1 should be consulted for a detailed dis- 
cussion of the use and limitations of this type of analysis. 
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It may also be worth, while' to study reference 2 "because, 
as early as 1886, Ayrton and Perry recognized the relation- 
ships that exist "between load and deflection of initially 
curved members. 



THEORY 



Consider the simple strut shown in figure 1. Assume 
the strut to be initially curved. Then under a load Pi, 
which is less than the critical load, it will have deflec- 
tions y x . These deflections can be accurately represented 
by the series 



TTX , . 2ttx . 
y"i = aj sxn — + ag sin + 

L L 



n = oo 



■I 



a n s in 



(1) 



n = i 



Now under an axial load P, which is greater than P x 
but less than the critical load, the deflections yj, will 
have been increased by amounts y a « If the total deflec- 
tions from the straight form for axial load P are y, 
then 



y = y z + y a 



(2) 



and the bending moment at any cross section is 



M = P(yj, + y a ) 



(3) 



Let the bending moment at any cross section be 
when the axial load is P 1 . Then 

M i = Pi7i 



(4) 



If M a is the increase in the bending moment as the axial 
load increases from P, to P» then 
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M 2 = M - M x 

= P(7i + 7 a ) - P x yi (5) 

or 

M a = (P - P x ) y x + Py a (6) 

The deflections y 3 caused by axial loads in excess of 
P x are determined from the differential equation 

EI --|L= . Ma - (?) 

dx 2 

Substitution of equation (l) in (6) and equation (6) in 
(7) gives, after division by EI, 



a n = 00 

d y 3 p / P P i \ V nirx . . 



dx EI EI EI ' / L 

n = l 



The general solution of this equation is 



y_ = A sin — + B cos - 
3 3 



n = oo 

n-nx 



+ a (l - *i) Y~-^- S in iffi (9) 

v P y £ n - a 

n = l 



where 



j = /— do) 
P 



P - 
a = (11) 

P crit 
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font - ^ (12) 

L 

In order to satisfy the end conditions (y 2 = 0, for 
x = 0 and for x = L) for any value of j , it follows 
that the constants of integration, A and B, must each 
equal zero. Therefore, the deflections y 2 are given by 
the equation 

n = 00 „ 

y = y - y = a U )> ~ sin 

3 i \ P J£ n - a L 

n = i 



■I 



n-rrx 

n_=_« ■ a n sin -j— 



n~=~i rS P crit " P * 1 



P - P, 



a, sin OS a s sin 

1 L L 



. + — (i 3 ) 

'crit ~ p i , 2 p crit " P i 



_1 ^ 



P - P p „ p 

As P approaches ^crit' f* rs * term in the series of 

equation (13) predominates. In this case it is possible 
to write as an approximation for equation (13) 

TTX 

a sin — 

1 L 

y 2 = y - y x = : U4) 

p crit - F i i 

p - p, 



But equation (14) can also be written in the form 
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v _ v a., i sin — — 

y - y, 7 - 7% v 



+ £- (15) 



? " P i P crit " P i P crit ' P i 



It is recalled tliat y - y x is the amount "by which 
the deflections are increased when the axial load on the 
strut is increased from P x to P. For any assumed ini- 
tial load r 1 the difference P cr it ~ ? i is a constant. 
Also, for any assumed cross section at which the lateral 

TTX 



a, sin -- 
L 



deflections 7 - 7 1 are measured, the term 

y _ y 

is a constant. Hence, if i- is plotted as ordinate 

against 7 - 7 X a s abscissa (fig. 2), it is recognized 
that equation (15) is a straight line. This line will cut 

the horizontal axis at the distance j^a x sin.~-J from the 



origin and the inverse slope of the line is P C rit " P i 4 

Thus if simultaneous readings of axial load and deflection 
are taken during a column test "beginning with any load Pi 
as the initial reading and these data are plotted in the 
manner just described, the reciprocal of the slope of the 
straight line obtained is the value of P C rit " p i • Tiie 
value of Pcrit ^ s tl3 - en .obtained from the relation 

p crit = ( p crit " p i) + p i < 16 > 

As mentioned by Southwell in reference 1, the main in- 
terest of this method of analysis lies in its generality 
because, in all ordinary examples of elastic instability, 
the same type of differential equation governs the deflec- 
tion as controlled by its' initial value, provided that both 
deflections are small. 

APPLICATION TO EXPERIMENTAL RESULTS 



In reference 1, Southwell applied his method of anal- 
ysis to the results of eight column tests made by T. von 
Karman ■ and published in 1909. In order to show that the 
more general equations of the present paper apply equally 
as well, these same data are reanalyzed in table I and fig 
ure 3. The method of least squares was used to establish 
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the best-f i tting straight line for each set of data plotted 
in figure S. This procedure was used in order that the 
personal equation would he eliminated in the manner fol- 
lowed by Southwell. 

In table II the results of— Southwell 1 s analysis and 
the analysis made in this report are compared. Inspection 
of the last two columns shows that the analysis made in 
this report predicts the critical load as closely as the 
analysis made by Southwell and that both predictions aro in 
close agreement with theory. 

The series of equation (l) gives the deflection curve 
under load P 1 . This series also gives the deflection 
curve when P x = 0 except that each coefficient a n has a 
value different than when P x > 0. If b n is substituted 
for a n when P x = 0, the relation between a n and b n 

is given by equation (10) of reference 1, which is, in the 
notation of this paper, 

a n = : — (17) 

1 1 

n3 *crit 



From this relation it is concluded that, as T 1 approaches 
^crit' first term in the series of equation (l) pre- 

dominates. Thus a x is a close approximation of the do- 
flection y x at the middle of - the strut. 

If the deflections y recorded by von Karman at the 
middle of the strut (see table I) represent the true deflec- 
tions from the condition of zero load, the value of a.,_ 

deduced from the bost-fitting line in figure 3 should bo in 
close agreement with the measured v aluo of y x . Inspection 
of table II shows that those values of a x and y 1 are 
not always in close agreement. The disagreement is not con- 
fined to the short columns but is also present in one of 
the long columns (strut 2). This fact rules out yielding 
of the material at high stresses as a possible explanation 
of the disagreement. 

One explanation of the agreement of a x and y x in 
table II in some cases and disagreement in other cases is 
as follows: When a 1 and y 1 agree, the strut was very 



N.A.C.A. technical Not© No. 658 



7 



nearly straight or the initial reading was taken at a very 
low load. When a r and y x disagreed, the strut had in- 
itial deflection and the initial reading was taken at other 
than zero load. 

This explanation was reached as a result of the fol- 
lowing reasoning. _ There must be some load on the column 
to hold it in the testing machine. Consequently, if the 
initial readings are taken at a load greater than zero, 
the deflections recorded are smaller than they would have 
been had the initial reading been taken at zero load. 
Thus, when the value of a x deduced from the best-fitting 
line in figure 3 is in disagreement with the value of y x 
recorded at the middle of the strut, y x ' should always be 
less than a x . This conclusion is verified by the values 
of a x and given in table II except for strut 4b. tn 

this strut a number of the readings at the lower loads are 
known to be out of line with the rest of the readings. 
(See fig. 8 of reference 1.) 

Another explanation of the disagreement between a x 
and y x in table II is as follows: Slight variations in 
the cross-sectional area and the material properties are 
possible in any strut. The effect of these variations ap- 
pears in the values of a x and P C rit deduced from the 

best-fitting line in figure 3. On a percentage basis, ai 
is much more sensitive than P Cr ^^. to variations of the 
type mentioned. 

If b x is the value of a x when Pj_ = 0, the values 

of b x that correspond to a^ are obtained from equation 
(17) with n = 1. These values of bj_ are listed in table 
II where .the values of b a deduced by Southwell are also 
tabulated. Inspection of the values of b x (N". A. C . A. ) re- 
veals that the largest initial deflections are found in 
those struts for which a x and y x disagree. This fact 
adds weight to the first explanation of the disagreement 
between, a. x and y 1# 

When Southwell estimated the critical load for the 
struts tested by von Karman, the slopes of the straight 
lines in figure 8 of reference 1 were determined by experi- 
mental points taken near the critical load. According to 
theory, the deflection approaches infinity when P ap- 
proaches ^crit* Although Southwell intended to plot y/P 

against y it may be that (y ± A ) /P was plotted against 



a 



N.A.C.A. Technical Note No, 658 



y ±- A » where A is a constant error in the measurement 
of y. When A is small and F is near p cr it» tJie 8S " 

timated value of ? cr i£ is very nearly the same in each 

case. It is therefore to be expected that Southwell's es- 
timate of the critical load should be as good as the esti- 
mate made in this report. (See table II.) 



CONCLUSIONS 



1. For the analysis of experimental observations in 
problems of elastic stability, it is found that the follow- 
ing equation holdsi 

ttx 

a sin - — 

7 - yi y - Yi 1 l 



P - p x p crit - P i p crit " P ! 



where 



P and y are the load and the corresponding deflec- 
tion, respectively. 

P x and y 1 are initial values of P and y, respec- 
tively. 

p crit is *ke critical value of P. 

a, sin -r- is a constant related to y. . 

Ii 1 

y — y 

2. The straight line obtained by plotting — as 

P - P, 

ordinate against y - y 1 as abscissa cuts the horizontal 

TTX I 

&x sin — from the origin and the 
L J 



axis at the distance 
inverse slope of the line is P C rit ~ P i 



3. For the experimental data examined, the critical 
load obtained from the slope of the straight line estab- 
lished by the test data was found to agr-ee well with the 
theoretical critical load. The values of , "however, 
did not always agree well with the value of y x observed 
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at the middle of the strut, indicating that the initial 
reading may not have "been taken at zero load, 

4. It is not always practicable to obtain the ini- 
tial reading of load and deflection at zero load. For 
this reason the method described herein for 'the analysis 
of experimental observations in problems of elastic sta- 
bility is more useful than methods previously described 
in which the initial reading must be at zero load. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , July 20, 1938. 
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APPENDIX 

Derivation of Equation (13) from the Equations 

Given in Reference 1 

If it is assumed that the deflections y Q at zero 
load are given "by the equation 



n = <» 

y 



■"v " nm 



o 

n = i 

then, according to reference 1, the deflections y at load 
P ar e 

JL2-SP 

y . N 2 sin (19) 



n = i 1 - 



The deflections y x at load P 2 are therefore 



n = co 

b 



, = > !» sin 212 (20 ) 

A / p L 



P, 

n = i i _ 1 



n a P 



Subtraction of equation (20) from equation (19) gives 
n = eo 



n = i 



! L_ !_ 



rf P crit n3 p crit 



nrrx 

h n sin 



n = e o b n sin SEi 

"vT 1 L 

= > (21) 

n = i n s Pcrit - p i p i 

a»MB«M^W^W^M^»aM i-| aM ^ — * ' ir ""' L r 11 



P " P i p crit 
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But 

■b 



a 2 (17) 

1 - ~li 

n * p crit 



With this gu"b st i tut ion, equation (21) agrees with equation 
(13) 

. nnx 

n = » a„ sin 

V n L 
y - y x = \ : : (13) 

t- t*i "D TD 



n = i n P crit " P i 



P - P x 



- 1 



Deflection at One Load Expressed in Terms 

of Deflection at Another Load 

If equation (l) gives the deflections at load Pj, , 
equation (13) givos the increase in deflections that result 
when the load is increased from P L to P. Consequently, 
addition of equations (1) and (13) gives for y, the de- 
flections at load P, 

y = 7 X + y 2 

n = °° a„ sin KTTX 



-I 



L 



n = i P - P x 
1 _ 1 



113 p crit - p x 

. rra: . Stt^ 
ax sin — a a sin 

L L 

+ -{-.... (22) 



P - P. P - P, 
1 1 i _ 1 

p crit " p i 2 p crit " P x 



As P approaches Pcrit* tlie f ir st term in the series 
of equation (22) predominates. Thus when P approaches 



12 



N.A.C.A. Technical Note No. 658 



P cr i_t * it is possible to write as an approximation for 
equation (22) 

a x sin 

y = (23) 

P - P. 



1 - 



P crit " p : 



If y x is accurately given by the first term of the series 
(1), y a and y are accurately given by the first terms 
of the series (13) and (22), respectively. Consequently, 
as long as the deflection curve remains a sine curve, the 
relations given by equations (14), (15), and (23) are exact, 
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Tablt 1 



Tabla I 

T. von Karman'e Struts 
[Hlld Steel: E » 2,170,000 kg/cm"] 



Strut 



P 

axial 
load 



(kg) 



deflection 
at nlddle 
of strut 1 

(in) 



P - Pi 



(kg) 



y - t x 

(am) 



P - P, 
mm per kg 



2,260 
3,020 
3,170 
3,320 
3,470 
3,620 



0.01 
.025 
.04 
.06 
.09 
.25 



760 
910 
1,060 
1,210 
1,360 



0.015 
.03 
.05 
.03 
.2* 



0.1974 x 10" 

.3297 

-4717 

.6612 
1.765 



4,520 
4,830 
5,130 
5,280 
5,430 



0.02 
.05 
.11 
.24 
.86 



310 
610 
760 
910 



O.03 
.09 
.22 
.84 



0.9677 
1.475 
2.895 
9.231 



3a 



3b 



6,030 
7,540 
8,290 
8,520 



•7,840 
8,140 
8,290 
8,445 
8,600 



0.01 
.03 
.11 
• 52 



1,510 
2.260 
2,490 



0.02 
.10 
.51 



0.02 
.05 
.07 
.11 
.21 



150 
305 
460 



0.02 
.06 
.16 



0.1325 
.4425 
2.048 



1.333 
1.967 
3-478 



•9,050 
•9,660 
10,260 
10,560 
10,710 
10,860 
11,010 
11,160 



0.02 
.025 
.03 
.07 
.10 

• 13 
.25 

• 73 



300 
450 
600 
750 
900 



0.04 
.07 
.10 
.22 
-70 



1.333 
1.556 
1.667 
2.933 
7.778 



4b 



♦3,020 
•4,530 
•6,030 
•7,540 
•8,300 
9,050 
9,805 
9,960 
10,110 
10,260 
10,410 
10,560 
10,710 
10,860 



0.03 
.05 
.07 
.09 
.12 
.15 
.23 
.26 
.29 

:8 
.52 

.71 
1.46 



755 
910 
1,060 
1,210 
1,360 
I.510 
1,660 
1,810 



0.08 
.11 
.14 
.18 
.26 

•1 
.56 

1.31 



1.060 
1.209 
1.321 
1.488 
1.912 
2.45a 
3-373 
7.238 



•9,050 
•10,560 
10,860 
11,160 
11,470 
11,770 
12,070 
12,370 
12,520 



0.01 
.03 
.05 
.07 
.10 
.15 
.22 
.50 
.45 



300 
610 
910 
1,210 
1,510 
1,660 



0.02 
.05 
.10 
.17 

•f 5 
.40 



0.6667 
.8197 
1.099 
1.405 
1.656 
2.410 



•10,560 
•12,070 
12,370 
12,670 
12,970 
13,270 
13,430 
13,580 



0.01 
.04 
.06 
.10 
.15 
.25 
.34 
.74 



300 
600 
900 
1,060 
1,210 



0.04 
.09. 
.19 
.28 
.68 



1.333 
1.50O 
2.111 
2.642 
5.620 



The data for these loads were rejeoted by Southwell on grounds that 
are stated at the beginning of paragraph 10 In reference 1. 
Consequently , the lowest load not so rejected la here taken as P . 



Table n 

Summary of Analyses Made of von Karman's TestB 



Strut 


Si 

recorded 
In test 

(mm) 


"i 

deduced 
from best- 
fitting 
line in 
figure 3 

(mm) 


deduced value 


P crit ( 8atimated ) 


P , .given 
crit 

by theoret- 
ical- 
formula 

(kg) 


P crit 
^crit 


(estimated) 
(theoretical) 


South- 
well 

(mm) 


N.A-C.A. 

(mm) 


Southwell 
(kg) 


N.A.C.A. 

(kg) 


Southwell 


N.A.O.A. 


1 


0.01 


0-0164 


0.005 


0.006 


3,712 


3,711 


3,790 


0.980 




0.980 


3 


.02 


.0602 


.005 


.011 


5,453 


5,495 


6,475 


.995 




1.004 


3a 


.01 


.0135 


.005 


.004 


8,590 


8,587 


8,645 


.994 




.993 


3b 


.05 


.0679 


.005 


.005 


8,758 


8,794 


B.610 


1.017 




1.020 


4a 


.03 


.0831 


.003 


.007 


11,230 


11,269 


10,980 


1,022 




1.025 


4o 


.15 


.1217 


.030 


.022 


11,090 


11,037 


10,930 


1.015 




1.011 


5 


.05 


.1350 


.010 


.023 


12,816 


13,095 


3J3,7BO 


1.003 




1.033 


6 ' 


.06 


.1304 


.010 


.014 


13,750 


13,833 


13,980 


.984 




.990 
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Figs. 1,2 



P Pi 




y 

Figure 1.- Baler strut. 




Figure 2.- Graph of equation (15) . 



